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ABSTRACT 

We have undertaken an unprecedentedly large 1.1 millimeter continuum survey of three nearby star 
forming clouds using Bolocam at the Caltech Submillimeter Observatory. We mapped the largest 
areas in each cloud at millimeter or submillimeter wavelengths to date: 7.5 deg^ in Perseus (Paper 
I), 10.8 deg^ in Ophiuchus (Paper II), and 1.5 deg^ in Serpens with a resolution of 31", detecting 
122, 44, and 35 cores, respectively. Here we report on results of the Serpens survey and compare 
the three clouds. Average measured angular core sizes and their dependence on resolution suggest 
that many of the observed sources are consistent with power-law density profiles. Tests of the effects 
of cloud distance reveal that linear resolution strongly affects measured source sizes and densities, 
but not the shape of the mass distribution. Core mass distribution slopes in Perseus and Ophiuchus 
(a = 2.1 ± 0.1 and a = 2.1 ± 0.3) are consistent with recent measurements of the stellar IMF, whereas 
the Serpens distribution is flatter (a = 1.6 ± 0.2). We also compare the relative mass distribution 
shapes to predictions from turbulent fragmentation simulations. Dense cores constitute less than 10% 
of the total cloud mass in all three clouds, consistent with other measurements of low star-formation 
efficiencies. Furthermore, most cores are found at high column densities; more than 75% of 1.1 mm 
cores are associated with Ay ^ 8 mag in Perseus, 15 mag in Serpens, and 20 — 23 mag in Ophiuchus. 
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1. INTRODUCTION 

Large-scale physical conditions in molecular clouds in- 
fluence the outcome of local star formation, including 
the stellar initial mass function, star-formation efficiency, 
and the spa tial distribution of stars within clouds (e.g. 
lEvanslll999( l. The physical processes that provide sup- 
port of molecular clouds and control the fragmentation 
of cloud material into star-forming cores remain a matter 
of debate. In the classical picture magnetic fields provide 
support and co llapse occurs via ambipolar diffusion (e.g. 
iShu et allHOTl . but many simulations now suggest that 
turbulence domi nates both support and f ragmentation 
(for a review see iMac Low fc Klessenl[200l . 

Dense prestellar and protostellar condensations, 
or cores (for definiti ons and an overview see 
iDi Francesco et all I2005D , provide a crucial link be- 
tween the global processes that control star formation 
on large scales and the properties of newly-formed stars. 
The mass and spatial distributions of such cores retain 
imprints of the fragmentation process, prior to signifi- 
cant influence from later protostellar stages such as mass 
ejection in outflows, core dissipation, and dynamical 
interactions. These cold (10 K), dense (n > 10^ cm~'^) 
cores are most easily observed at millimeter and sub- 
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millimeter wavelengths where continuum emission from 
cold dust becomes optically thin and traces the total 
mass. Thus, complete maps of molecular clouds at 
millimeter wavelengths are important for addressing 
some of the outstanding questions in star formation. 

Recent advances in millimeter and submillimeter wave- 
length continuum detectors have enabled a number 
of large-scale surveys of nearby molecular clouds (e.g 
Johnstone et all 120041: iKirk et al . 2005; HatchcU et al 



2005: 'Eno ch et al.ir200a IStanke et al.ii2006i: Young et al 
2006). In addition to tracing the current and fu- 
ture star-forming activity of the clouds on large scales, 
millimeter and submillimeter observations are essential 
to understanding the properties of starless cores and 
the envelopes of the most deeply embedded protostars 
(for more on the utility of millimeter observations, see 
lEnoch et al.1l2006D . 

We have recently completed 1 . 1 mm surveys of Perseus 



( Enoch et al 



(l Young et al 



l2006l hereafter Paper I) and Ophiuchus 
|2006l hereafter Paper II) with Bolocam at 
the Caltech Submillimeter Observatory (CSO). In this 
work we present a similar 1.1 mm survey of Serpens, com- 
pleting our three-cloud study of nearby northern star- 
forming molecular clouds. Unlike previous work, our sur- 
veys not only cover the largest area in each cloud to date, 
but the uniform instrumental properties allow a compre- 
hensive comparison of the cloud environments in these 
three regions. A comparison of the results for all three 
clouds provides insights into global cloud conditions and 
highlights the influence that cloud environment has on 
properties of star forming cores. 

Background facts on the Perseus and Ophiuchus molec- 
ular clouds are discussed in Papers I and II. The Serpens 
molecular cloud is an active star formation region at a 
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distance of d = 260 ± 10 pc (jStraizys et al.lll996f ). Al- 
though the cloud exte nds more then 10 deg^ as mapped 
by optical extinction (jCambresvl I1999D . most observa- 
tions of the region have been focused near the main Ser- 
pens cluster at a Right Ascension (R.A.) of IS^'SO™ and 
declination (decl.) of 1°15' (J2000). 

The Serpens cluster is a highly extincted region with 
a high density of young stellar objects (YSOs), includ- 
ing a number of Class protostars. It has been studied 
extensively at near-infrared, far-inf rared, subm illimcter, 
and millimeter wavelengths (e.g. lEiroa fc Casali ,1992i: 
Hurt fc Barsonvl 119961: iLarsson et all 120001: iDavis et al l 



1999| iCasah. Eiroa. fc DuncanI 119931 iTesti fc SargartI 
19981). Some recent work has also drawn attention to a 



less well known cluste r to the south, sometimes r eferred 
to as Serpens/G3-G6 (|Diupvik et al.| [2006: Harvev et al.l 
|2006[) . Beyond these two clusters no continuum millime- 
ter or submillimcter continuum surveys have been done 
that could shed light on large-scale star formation pro- 
cesses. 

Following Papers I and II, we utilize the wide-field 
mapping capabilities of Bolocam, a large format bolome- 
ter array at the CSO, to complete millimeter continuum 
observations of 1.5 deg^ of the Serpens cloud. These ob- 
servations are coordinated to cover the area mapped with 
Spitzer Space Telescope IRAC and MIPS observations 
of Se rpens from the "Cores to Disks" (c2d I Evans et al.l 
|2003() legacy project. While millimeter and submillimc- 
ter observations are essential to understanding the prop- 
erties of dense prestellar cores and protostellar envelopes, 
infrared observations are necessary to characterize the 
protostars embedded within those envelopes. In a future 
paper (M. Enoch et al. 2007, in preparation) we will take 
advantage of the overlap between our 1.1 mm maps and 
the c2d Spitzer Legacy maps to characterize the deeply 
embedded and prestellar populations in Serpens. 

Here the Bolocam observations and data reduction 
(§ HI? general cloud morphology and source properties 
(§ El, and summary (§ |4]) for Serpens are briefly pre- 
sented. In § [5] we compare the millimeter survey results 
for Serpens, Perseus, and Ophiuchus. First, we outline 
our operational definition of a millimeter core including 
instrumental effects in § 15.11 and discuss the observa- 
tional biases introduced by different cloud distances in 
§ 15.21 Physical implications of source sizes and shapes, 
and differences between the clouds, are discussed in 15.31 
We examine source densities and compare the cloud mass 
versus size distributions in § 15.41 We analyze the core 
mass distributions and their relation to cloud turbulence 
in § 15.51 and the spatial distributions of the core samples 



and clustering in § 15.61 The relationship of dense cores 
to the surrounding cloud column density and the mass 
fraction in dense cores are discussed in § 15.71 and § 15.81 
We end with a summary in § [6l 

2. SERPENS OBSERVATIONS AND DATA REDUCTION 

Observations and data reduction for Serpens follow the 
same methodology as for Perseus (Paper I) and Ophi- 
uchus (Paper II). The data reduction techniques we have 
developed for the molecular cloud data are described in 
detail in Paper I and Paper II, including removal of sky 
noise, construction of pointing and calibration models, 
application of iterative mapping, and source extraction. 
Given these previous descriptions, only details specific 



to Serpens will be presented here. Further information 
about the Boloc am instrument and r eduction pipeline is 
also available in lLaurent et al.l ()2005l ). 

2.1. Observations 

As for Perseus and Ophiuchus, millimeter continuum 
observations of Serpens were made with Bolocam^ at the 
CSO on Mauna Kea, Hawaii. Bolocam is a 144-element 
bolometer array that operates at mil limeter wavel engths 
and is designed to map large fields (| Glenn et al .1 120031. 
Observations of Serpens were carried out at a wavelength 
of A = 1.1 mm, where the field of view is 7.'5 and the 
beam is well approximated by a Gaussian with full- width 
at half-maximum (FWHM) size of 31". The instrument 
has a bandwidth of 45 GHz at A = l.lmni, excluding the 
CO (J=2-l) line to approximately 99%. 

The 1.1mm observations were designed to cover a re- 
gion with Ay > 6 magnitudes (Ay > 6™) in the visual 
extinction map of ICambresvl (|1999( ) , shown in Figure [TJ 
As demonstrated in Figure [1] this ensures that the Bolo- 
cam observations overlap as closely as possible with the 
region of Serpens observed with Spitzer IRAC and MIPS 
as part of the c2d Legacy project. In practice, the Bolo- 
cam survey covers a region slightly larger than the IRAC 
map, and slightly smaller than the MIPS map. 

Serpens was observed during two separate runs in 2003 
May 21-June 9 and 2005 June 26-30. During the 2003 
run, 94 of the 144 channels were operational, compared 
to 109 during the 2005 run. Scans of Serpens were made 
at a rate of 60" sec~^ with no chopping of the secondary. 
The final map consists of 13 scans from 2003 in good 
weather, and 17 scans from 2005 in somewhat poorer 
conditions. Each scan covered the entire 1.5 deg^ area 
and took 35 — 40 minutes to complete depending on the 
scan direction. Scans were made in two orthogonal di- 
rections, approximately half in R.A. and half in decl. 
This strategy allows for good cross-linking in the final 
map, sub-Nyquist sampling, and minimal striping from 
1// noise. Scans were observed in sets of three offset by 
—43", 0", -1-43" to optimize coverage. 

In addition, small maps of pointing sources were ob- 
served approximately every 2 hours, and at least one pri- 
mary flux calibration source, including Neptune, Uranus, 
and Mars, was observed each night. Several larger beam 
maps of planets were also made during each run, to char- 
acterize the Bolocam beam at 1.1 mm. 

2.2. Pointing and Flux Calibration 

A pointing model for Serpens was generated using 
two nearby pointing sources, G 34.3 and the quasar 
1749-f096. After application of the pointing model, a 
comparison to the literature S CUBA 850 /xm p ositions 
of four bright known sources (IDavis et al.l 11993 ) in the 
main Serpens cluster indicated a constant positional off- 
set of ((5R.A., Mecl.) = (5", -10"). We corrected for the 
positional offset, but estimate an uncertainty in the ab- 
solute pointing of 10", still small compared to the beam 
size of 31". The relative pointing errors, which cause 
blurring of sources and an increase in the effective beam 
size, should be much smaller, approximately 5". Relative 
pointing errors are characterized by the rms pointing un- 

^ jhttp: //www. CSO ■ caltech.edu/bolocani] 
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certainty, derived from the deviations of G 34.3 from the 
pointing model. 

The Bolocam flux cahbration method makes real-time 
corrections for the atmospheric attenuation and bolome- 
ter operating point using the bolometer optical loading, 
by calculating the calibration factor as a function of the 
bolometer DC resistence. Calibrator maps of Neptune, 
Uranus, and G 34.3, observed at least once per night, 
were used to construct a calibration curve for each run. 
A systematic uncertainty of approximately 10% is as- 
sociated with the absolute flux calibration, but relative 
fluxes should be much more accurate. 

2.3. Cleaning and Iterative mapping 

Aggressive sky subtraction techniques are required for 
Bolocam data to remove sky noise, which dominates 
over the astronomical signal before cleaning. As in Pa- 
pers I and II, we remove sky noise from the Serpens 
scan s using Princip a l Com ponent Analysis (PCA) clean- 
ing (jLaurent et al.l ()2005l ) and references therein), by 
subtracting 3 PCA components. 

As described in Paper I, PCA cleaning removes some 
source flux from the map as well as sky noise, necessitat- 
ing the use of an iterative mapping procedure to recover 
lost astronomical flux. In this procedure, sky subtraction 
is refined by iteratively removing a source model (derived 
from the cleaned map) from the raw data, thereby reduc- 
ing contamination of the sky template by bright sources. 
Simulations show that at least 98% of source flux density 
is recovered after 5 iterations, with the exception of very 
large (> 100" FWHM) faint sources, for which we only 
recover 90% or less of the true flux density (See Paper I) . 
We conservatively estimate a 10% residual photometric 
uncertainty from sky subtraction after iterative mapping, 
in addition to the 10% uncertainty in the absolute flux 
calibration. 

Data from the 2003 and 2005 observing runs were iter- 
atively mapped separately because they required differ- 
ent pointing and flux calibration models. After iterative 
mapping the two epochs were averaged, weighted by the 
square root of the observational coverage. 

2.4. Source Identification 

Millimeter cores are identified in an optimally filtered 
map using the source extraction method described in Pa- 
per I. Because the signal from a point source lies in a 
limited frequency band, we can use an optimal (Wiener) 
filter in Fourier space to attenuate 1/f noise at low fre- 
quencies, as well as high frequency noise above the sig- 
nal band. The optimal filter preserves the resolution of 
the map and the peak brightness of point sources but 
reduces the rms noise per pixel by approximately -v/S, 
thereby optimizing the source signal-to- noise (S/N). Ex- 
tended sources will have slightly enhanced peak values in 
the optimally filtered map. After optimal filtering, the 
map is trimmed to remove areas of low coverage. Note 
that the optimally filtered map is used for source detec- 
tion only; all photometry is measured in the unfiltered 
map, and all maps displayed here are unfiltered. 

Observational coverage, which depends on the scan 
strategy, number of scans, and number of bolometer 
channels, was very uniform for Serpens; trimming re- 
gions where the coverage was less than 30% of the peak 



coverage was equivalent to cutting off the noisy outer 
edges of the map. The average coverage of the Serpens 
map is 1600 hits per pixel, where a hit means a bolome- 
ter passed over this position, with differences in hits per 
pixel across the map of 18%. The average coverage cor- 
responds to an integration time of 13 minutes per pixel, 
although individual pixels are not independent because 
the map is over-sampled. 

For each pixel, the local rms noise is calculated in small 
(45 arcmin^) boxes using a noise map from which sources 
have been removed. This noise map is derived as part 
of the iterative mapping process (see Paper I). A sim- 
ple peak-finding routine identifies all pixels in the op- 
timally filtered map more than 5a above the local rms 
noise level, and source positions are determined using 
an IDL centroiding routine. Each new source must be 
separated by at least a beam size from any previously 
identified source centroid. The centroid is a weighted 
average position based on the surface brightness within 
a specified aperture, and is computed as the position 
at which the derivatives of the partial sums of the in- 
put image over (y,x) with respect to (x,y) equal zero. A 
given centroid is considered "well defined" as long as the 
computed derivatives are decreasing. All sources were 
additionally inspected by eye to remove spurious peaks 
near the noisier edges of the map. 

3. SERPENS RESULTS 

The final 10" pixel""'^ Serpens map is shown in Figure[2l 
with the we ll known northern Serpens cluster. Cluster A 
(|Harvev et al. 2006), indicated, as well as the southern 
cluster. Cluster B (Serpens/G3-G6). Covering a total 
area of 1.5 deg^, or 30.9 pc^ at a distance of 260 pc, the 
map has a linear resolution of 7.3 x 10"^ AU. 

We identify 35 sources above the 5cr detection limit in 
the Serpens map. The 5cr limit is based on the local rms 
noise, and is typically 50 mJy beam~^. The noise map 
for Serpens is shown in Figure [31 with the positions of 
identified sources overlaid. The mean rms noise is 9.5 
mJy beam"""^, but is higher near bright sources. Most 
of the 18% variations in the local rms noise occur in the 
main cluster region, where calculation of the noise is con- 
fused by residual artifacts from bright sources. Such ar- 
tifacts must contribute significantly to noise fluctuations; 
coverage variations alone would predict rms noise varia- 
tions of only \/T8 — 4%. This means that faint sources 
near bright regions have a slightly lower chance of being 
detected than those in isolation. 

Source positions are listed in Table [1] and identified 
by red circles in Figure ID Figure [4] also shows magni- 
fications of the more densely populated source regions, 
including Cluster A and Cluster B. We do not see any 
circularly symmetric extended emission on scales > 3' in 
the map. It should, in principle, be possible to recover 
symmetric structures up to the array size of 7f5, but our 
simulations show that sources > 4' in size are severely af- 
fected by cleaning and therefore difficult to fully recover 
with iterative mapping. The map does contain larger fil- 
amentary structures up to 8' long. In particular, the long 
filament between Cluster A and Cluster B is reminiscent 
of the elongated ridge near Bl in Perseus (Paper I). The 
Serpens filament does not contain the bright compact 
sources at either end that are seen in the Bl ridge, how- 
ever. 
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Previous millimeter-wavelength maps of Clus- 

ter A, such as the lA mm UKT14 map of 

ICasali. Eiroa. &: DuncanI ( 19931) and the 850 fim 
SCUBA map of lDavis et al ( 1999D . generally agree with 
our results in terms of morphology and source structure. 
Not all of the individual 850 /xm sources are detected 
by our peak finding routine, presumably due to the 
poorer resolution of Bolocam (30") compared to SCUBA 
(14"), but most can, in fact, be identified by eye in the 
Bolocam map. An IRAM 1.3 mm continuum map of 
Cluster B with 11" resolution (jPiupvik et all |2006l ) is 
visually quite similar to our Bolocam map of the region. 
We detect each of the four 1.3 mm sources identified 
(MMSl-4), although the 1.3 mm triplet MMSl is seen 
as a single extended source in our map. 

Most of the brightest cores, in particular those in 
Cluster A, are associated with known YSOs inc l udin ; 
a number of Class objects (|Hurt fc Barsonvl 1199 
iHarvev et all Il006). All bright 1.1 mm sources are 
aligned with bright 160 fim emission in the Spitzer MIPS 
map of Serpens o bserved by the c2d legacy project 
(jHarvev et al.|[2007l ). Fainter millimeter sources are usu- 
ally associated with extended 160 fim filaments, but do 
not necessarily correspond to point sources in the MIPS 
map. Conversely, one bright extended region of 160 /im 
emission just south of Cluster A contains no 1.1 mm 
sources. This area also exhibits extended emission at 70 
and 24 /im, which may be indicative of warmer, more 
diffuse material than that of the dense cores detected at 
1.1 mm. 

Despite the low rms noise level achieved in Serpens, 
very few sources are seen outside the main clusters; most 
of the area that we mapped appears devoid of 1.1 mm 
emission, despite being in a region of high extinction. 
Figure [5] shows a comparison between the Bolocam mil- 
limeter map (gray-scale) and visual extinction (contours) 
derived from c2d near- and mid-infrared Spitzer data. 

The majority of sources (90-95%) detected by IRAC 
and MIPS in the c2d clouds have spectral energy distri- 
butions characteristic of reddened stars. Thus we have 
measures of the visual extinction for many lines of sight 
through the molecular clouds imaged by c2d. Line-of- 
sight extincti on values are derived by fitting the Rv=5.5 
dust model of lWeingartner fc Draind (l2001h to the near - 
infrared through mid-infrared SED (jEvans "etal][200l . 
For each of the three clouds, the derived line-of-sight 
extinctions were convolved with uniformly spaced 90" 
Gaussian beams to construct an extinction map. 

Extinction maps used throughout this paper for Ser- 
pens, Perseus, and Ophiuchus are derived from c2d data 
by this method. The c2d extinction maps accurately 
trace column densities up to Ay ^ 40™, but are rela- 
tively insensitive to small regions of high volume den- 
sity, because they rely on the detection of background 
stars. Thus the extinction maps are complementary to 
the 1.1 mm observations, which trace high volume den- 
sity structures (see § 15. ip . In Figure EJ the Ay map for 
Serpens is smoothed to an effective resolution of 2.'5. 

As can be seen in Figure O nearly all Serpens mil- 
limeter sources lie within regions of high visual extinc- 
tion {Ay > 10™) and, in particular, all bright Bolocam 
sources are associated with areas of Ay > 15™. Nev- 
ertheless, there are a number of high extinction areas 
{Ay > 12™) with no detectable 1.1 mm sources. A simi- 



lar general trend was noted in both Perseus (Paper I) and 
Ophiuchus (Paper II), with relatively few sources found 
outside the major groups and clusters associated with 
the highest extinction. § 15.71 examines the relationship 
between Ay and 1.1 mm sources in more detail. 

3.1. Source Properties 
3.1.1. Positions and Photometry 

Positions, peak flux densities, and S/N for the 35 
1.1 mm sources identified in the Bolocam map of Ser- 
pens are listed in Table [T] The S/N ratio is measured 
in the optimally filtered map, whereas photometry and 
all other source properties are measured in the unfil- 
tered, surface brightness normalized map. The peak flux 
density per beam (7^) is given in mJy beam~^ (1 mJy 
beam~^ = 0.04 MJy sr~^). Uncertainties in photome- 
try are calculated from the local rms beam^-'^, calculated 
as in § 12.41 An additional systematic error of 15% is 
associated with all flux densities, from the absolute cal- 
ibration uncertainty and the systematic bias remaining 
after iterative mapping. Table [T] also lists the most com- 
monly used name from the literature for known sources, 
and indicates if the 1.1 mm source is coincident (within 
60") with a MIPS 2 4 /im source from the c2d database 
(jHarvev et al.ll2007t ). 

Table [2] lists photometry in flxed apertures of diameter 
40", 80", and 120", and the total integrated flux den- 
sity {S^). Integrated flux densities are measured assum- 
ing a sky value of zero, and include a correction for the 
Gaussian beam so that a point source has the same in- 
tegrated flux density in all apertures. No integrated flux 
density is given if the distance to the nearest neighbor- 
ing source is smaller than the aperture diameter. The 
total flux density is integrated in the largest aperture 
(30" — 120" diameters in steps of 10") that is smaller 
than the distance to the nearest neighboring source. Un- 
certainties are aap = <Jmb{Qapl()mb), whcrc Omb is the 
local rms beam~^ and (^op, Qmb) are the aperture and 
beam FWHM respectively. 

Peak and total flux density distributions for the 35 
1.1 mm sources in Serpens are shown in Figure [51 with 
the 5cr detection limit indicated. In general, source total 
flux densities are larger than peak flux densities because 
most sources in the map are extended, with sizes larger 
than the beam. Both distributions look bimodal; all of 
the sources in the brighter peak are in either Cluster A or 
Cluster B. The mean peak flux density of the sample is 
0.5 Jy beam~^, and the mean total flux density is 1.0 Jy, 
both with large standard deviations of order the mean 
value. Peak Ay values of the cores, calculated from the 
peak flux density as in § 13.1.31 are indicated on the upper 
axis. 

3.1.2. Sizes and Shapes 

Source FWHM sizes and position angles (PA, mea- 
sured east of north) are measured by fltting an elliptical 
Gaussian after masking out nearby sources using a mask 
radius equal to half the distance to the nearest neighbor. 
The best fit major and minor axis sizes and PAs are listed 
in Table [2) Errors given are the formal fitting errors and 
do not include uncertainties due to residual cleaning ef- 
fects, which are of order 10 - 15% for the FWHM and 5° 
for the PA (see Paper I). 
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As can be seen in Figure [7l the minor axis FWHM 
values are fairly narrowly distributed around the sample 
mean of 49", with a standard deviation of 12". The 
major axis FWHM have a similarly narrow distribution 
with a mean of 63" and a scatter of 12". On average 
sources are slightly elongated, with a mean axis ratio at 
the half-max contour (major axis FWHM / minor axis 
FWHM) of 1.3. 

A morphology keyword for each source is also given 
in Tabled! to describe the general source shape and en- 
vironment. Keywords indicate if the source is multiple 
(within 3' of another source), extended (major axis full- 
width at 2(7 > 1'), elongated (axis ratio at 4a > 1.2), 
round (axis ratio at 4<7 < 1.2), or weak (peak flux densi- 
ties less than 5 times the rms per pixel in the unfiltered 
map). The majority (28/35) of sources are multiple by 
this definition, and nearly all (32/35) are extended at the 
2(7 contour. 

3.1.3. Masses, Densities, and Extinctions 

The total mass M of gas and dust in a core is pro- 
portional to the total flux density S^, assuming the dust 
emission at 1.1 mm is optically thin and both the dust 
temperature and opacity are independent of position 
within a core: 

M = \ , (1) 

where Ki.imm — 0.0114 cm^ g^^ is the dust opacity per 
gram of gas, d = 260 pc is the distance, and To = 10 K 
is the dust temperature. A gas to dust mass ratio of 
100 is incl uded in KiAmm- The dust op acity is interpo- 
lated from lOssenkopf &: Hennin3 (|1994 Table 1 column 
5), for dust grains with thin ice mantles, coagulated for 
10^ years at a gas density of 10^ cm~^. This opacity has 
been found to b e the best fit in a number of radiative 
transfer models p vanseLaU |200l|; IShirlev et aP l2002t 
lYoung et al.ll2003f K 

Masses calculated in this manner, assuming a dust 
temperature of To = 10 K for all sources, are listed in 
Table [2j Uncertainties given are from the uncertainty in 
the total fiux density only; additional uncertainties from 
K, Td, and d together introduce a total uncertainty in the 
mass of up to a factor of 4 or more (also see Paper I) . The 
dust opacity is uncertain by up to a factor of two or more 
([Osscnkopf & Hcnning 1994), owing to large uncertain- 
ties in the assumed dust properties, and possible posi- 
tional variations of within a core. Variations in both 
Ki, and the cloud distance d have smaller effects than 
dust temperature uncertainties for the range of plausible 
values (ki irnm = 0.005-0.02 cm^ g-i; d = 200-300 pc; 
Ti, = 5 - 30 K). 

Assuming To = 10 K is a reasonable compromise to 
cover both prestellar and protostellar sources, based on 
the results of more de tailed radiativ e transfer models 
(|Evans et al. 2001; Shirl ev et all2002l see also discussion 
in Paper I). A temperature of 10 K will result in overes- 
timates of the masses of protostellar cores, which will be 
warmer on the inside, by up to a factor of three (for a 
20 K source) . Temperature errors are not as problematic 
as they might seem, however, as most of the envelope 
mass is located at large radii and low temperatures (see 
also Paper I). 

The total mass of the 35 1.1 mm cores in Serpens is 



92Mq, which is only 2.7% of the total cloud mass (3470 
Mq). The cloud mass is estimated from the c2d visual ex- 
tinction map using iV(H2)/Ay — 0.94 x 10^^ mag cm~^ 
(|Bohlin et al.| [l9781 and 

Af (cloud) = d^mHlJ.H^VL^N{Yl2) (2) 

where d is the distance, f2 is the solid angle, ^ indicates 
summation over all Ay > 2™ pixels in the extinction 
map, and fiH2 = 2.8 is the mean molecular weight per 
H2 molecule. 

Figure [8] shows the differential mass function of all 
1.1 mm sources in Serpens. The point source detec- 
tion limit of O.I3M0 is indicated, as well as the 50% 
completeness limit for sources of FWHM 55" (0.35Mq), 
which is the average size of the sample. Completeness is 
determined from Monte Carlo simulations of simulated 
sources inserted into the raw data and run through the 
reduction pipeline, as described in Paper I. The best fit 
power law slope to dN/dM oc M~" is shown {a = 1.6), 
as well as the best fit lognormal slope. 

Our ma ss distribution has a flatter slope than that 
found by Testi fc Sargentl ()1998f ) from higher resolution 
(5") OVRO observations (a = -2.1 for M > 0.35Mq). 
This may be due in part to the fact that most of our 
detections, at least 25/35, lie outsi de the 5f5 x 5.'5 area 
observed bv iTesti fc Sargentl (|1998f l. The resolution dif- 
ferences of the observations may also contribute signifi- 
cantly; for example, a number of our bright sources break 
down into multiple objects in the 5" resolution map. 

From the peak flux density Ii, we calculate the central 
H2 column density for each source: 

^(^2) = „ (3) 

Here r2„ifc is the beam solid angle, uih is the mass of hy- 
drogen, Ki.imm = 0.0114 cm^g~^ is the dust opacity per 
gram of gas, B^, is the Planck function, Td=10 K is the 
dust temperature, and = 2.8 as above. From column 
density we conv ert to ex tinction: Ay = A^(H2)/0.94 x 
10^1 mag cm-2 (jBohlin e t al. 197j), adopting Ry = 3.1. 
We note, however, that this relation was determined for 
the diffuse interstellar medium and may not be ideal for 
the highly extincted lines of sight probed here. 

The resulting peak Ay values are listed in Table [H 
The average peak Ay of the sample is 41™ with a large 
standard deviation of 55™ and a maximum Ay of 256™. 
Extinctions calculated from the millimeter emission are 
generally higher than those from the c2d visual extinction 
map by approximately a factor of 7, likely a combination 
of both the higher resolution of the Bolocam map (30" 
compared with 90") and the fact that the extinction map 
cannot trace the highest volume densities because it relies 
on the detection of background sources. Grain growth 
in dense cores beyond that included in the dust opacity 
from Osscnkopf & Henning (1994) could also lead to an 
overestimate of the Ay from out 1.1 mm data. 

Also listed in Table [2] is the mean particle density: 

where M is the total mass, R is the linear deconvolved 
half-width at half-max size, and fip — 2.33 is the mean 
molecular weight per particle. The median of the source 
mean densities is 2.3 x 10^ cm~'^, with values ranging 
from 3.1 X 10^ to 6.1 x 10^ cm~3. 
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4. SERPENS SUMMARY 

We have completed a 1.1 mm dust continuum survey 
of Serpens, covering 1.5 deg^, with Bolocam at the CSO. 
We identify 35 1.1 mm sources in Serpens above a 5a 
detection Hmit which is 50 mJy beam~^, or O.ISMq, on 
average. The sample has an average mass of 2.6M0, 
and an average source FWHM size of 55". On average, 
sources are slightly elongated with a mean axis ratio at 
half- max of 1.3. The differential mass distribution of all 
35 cores is consistent with a power law of slope a = 
1.6 ±0.2 above O.35M0. The total mass in dense 1.1 mm 
cores in Serpens is 92Mq, accounting for 2.7% of the total 
cloud mass, as estimated from our c2d visual extinction 
map. 

5. THREE-CLOUD COMPARISON OF PERSEUS, 
OPHIUCHUS, AND SERPENS 

The survey of Serpens completes a three-cloud 
study investigating the properties of millimeter emis- 
sion in nea rby star-forming mo lecular clouds: Pe rseus 
(jEnoch et a l. 2006), Ophiuchus (|Young et al.ll2006f ). and 
Serpens. Having presented the results for the Serpens 
cloud, we now compare the three clouds. 

Our large-scale millimeter surveys of Perseus, Ophi- 
uchus, and Serpens, completed with the same instrument 
and reduction techniques, provide us with a unique ba- 
sis for comparing the properties of 1.1 mm emission in 
a variety of star-forming environments. In the following 
sections we examine similarities and differences in the 
samples of star-forming cores, and discuss implications 
for physical properties of cores as well as global cloud 
conditions. 

5.1. What is a core? 

Before comparing results from the three clouds, we first 
describe our operational definition of a millimeter core. 
The response of Bolocam to extended emission, together 
with observed sensitivity limits determine the type of 
structure that is detectable in our 1.1 mm maps The 
Bolocam 1.1 mm observations presented here are sensi- 
tive to sub-structures in molecular clouds with volume 
density n > 2 x 10"* cm~^. One way to see this is to cal- 
culate the mean density along the curve defined by the 
detection as a function of size in each cloud. 

Figure [9] demonstrates how the completeness in each 
cloud varies as a function of source size. Plotted sym- 
bols give the total mass versus linear deconvolved size for 
all sources detected in each cloud, and lines indicate the 
empirically derived 50% completeness limits. For Ophi- 
uchus the average completeness curve is plotted; as the 
rms noise varies considerably in Ophiuchus, some regions 
have higher or lower completeness limits than the curve 
shown here. Completeness is determined from Monte 
Carlo simulations by adding simulated sources to the raw 
data, processing them in the same way as the real data, 
and attempting to detect them using our peak-finding 
algorithm (see Paper I). We are biased against detecting 
large diffuse sources because we detect sources based on 
their peak flux density, whereas the mass is calculated 
from the total flux, which scales approximately as the 
size squared. 

Calculating mean densities along the 50% complete- 
ness curve in each cloud yields num ~ 3 — 7 x lO'' cm^'^ 
in Serpens, num ~ 2 — 4 x 10** cm~^ in Perseus, and 



niim ^ 10 — 30 X 10^ cm~^ in Ophiuchus. By compari- 
son, the mean cloud density as probed by the extinction 
map is approximately 1000 cm~'^ in Serpens, 220 cm^'^ 
in Perseus, and 390 cm~'^ in Ophiuchus. To be identified 
as a core, therefore, individual structures must have a 
mean density n > 2 x 10"* cm~^, and a contrast com- 
pared to the average background density of at least 30- 
100. The mean cloud density is estimated from the to- 
tal cloud mass (§ I3.1.3P and assumes a cloud volume of 
V = A^-^, where A is the area of the extinction map 
within the Ay = 2 contour. 

Although we are primarily sensitive to cores with high 
density contrast compared to the background, it is clear 
that there is structure in the 1.1 mm map at lower con- 
trasts as well, and that many cores are embedded within 
lower density filaments. The total mass in each of the 
1.1 mm maps, calculated from the sum of all pixels 
> 5(7, is approximately twice the mass in dense cores: 
176Mo versus 92Mq in Serpens, 376M0 versus 278Mq 
in Perseus, and 83Mq versus 44M0 in Ophiuchus, for 
ratios of total 1.1 mm mass to total core mass of 1.9, 
1.4, and 1.9 respectively. Thus about half the mass de- 
tectable at 1.1 mm is not contained in dense cores, but is 
rather in the "foothills" between high density cores and 
the lower density cloud medium. 

Structures that meet the above sensitivity criteria and 
are identified by our peak-finding routine are considered 
cores. Our peak-finding method will cause extended fil- 
aments to be broken up into several separate "cores" if 
there are local maxima in the filament separated by more 
than one beam size, and if each has a well-defined cen- 
troid (see § 12. 4p . There is some question as to whether 
these objects should be considered separate sources or 
a single extended structure, but we believe that our 
method is more reliable for these d ata than alternative 
methods such as Clumpfind (Willia ms et al.i ri994). In 
Paper I we found that faint extended sources in our maps, 
which one would consider single if examining by eye, 
are often partitioned into multiple sources by Clumpfind. 
Using our method, one filamentary structure in Serpens 
is broken up into several sources, as are two filaments in 
Perseus. 

Monte Carlo tests were done to quantify biases and 
systematic errors introduced by the cleaning and itera- 
tive mapping process, which affect the kind of structure 
we can detect. Measured FWHM, axis ratios, and posi- 
tion angles are not significantly affected by either clean- 
ing or iterative mapping for sources with FWHM< 120"; 
likewise, any loss of flux for such sources has an am- 
plitude less than that of the rms noise. Sources with 
FWHM> 200" are detected, but with reduced flux den- 
sity (by up to 50%), and large errors in the measured 
FWHM sizes of up to a factor of two. 

The limitation on measurable core sizes of approxi- 
mately 120" corresponds to 3 x lO'' AU in Perseus and 
Serpens, and 1.5 x lO'' AU in Ophiuchus. Note that 
these sizes are of order the median core separation in 
each cloud (see § 15. 6|) . meaning that we are just as likely 
to be limited by the crowding of cores as by our sensi- 
tivity limits in the measurement of large cores. The de- 
pendence of measurable core size on cloud distance can 
be seen in Figure [9l where the completeness rises steeply 
at smaller linear deconvolved sizes for Ophiuchus than 
for Perseus or Serpens. Thus we are biased against mea- 
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suring large cores in Ophiuchus compared to the other 
two clouds. Although there are cores in our sample with 
sizes up to 3 X 10^ AU (Figure (HI), most cores have sizes 
substantially smaller than the largest measurable value, 
again indicating that we are not limited by systematics. 

To summarize, the Bolocam 1.1 mm observations pre- 
sented here naturally pick out sub-structures in molecu- 
lar clouds with high volume density (n > 2 x 10^ cm~^). 
These millimeter cores have a contrast of at least 30—100 
compared to the average cloud density as measured by 
the visual extinction map (2 — 10 x 10'^ cm~'^). Many 
cores are embedded in lower density extended structures, 
which contribute approximately half the mass measur- 
able in the Bolocam 1.1 mm maps. Finally, Monte Carlo 
tests indicate that we can detect cores with intrinsic sizes 
up to approximately 120". 

5.2. Distance effects 

To test the effects of instrumental resolution and its de- 
pendence on distance, we convolve the Ophiuchus map 
with a larger beam to simulate putting it at approxi- 
mately the same distance as Perseus and Serpens. After 
convolving the unfiltered Ophiuchus map to 62" reso- 
lution, we apply the optimal filter and re-compute the 
local rms noise, as described for Serpens (§ 12. 4[) . The 
pixel scale in the convolved map is still 10" pixel^^, but 
the resolution is now 62" and the rms is lower than in the 
original map, with a median value of 17 mJy beam~^ in 
the main L 1688 region. Source detection and photom- 
etry is carried out in the same way as for the original 
map, with the exception that a 62" beam is assumed. 

We detect 26 sources in the degraded-resolution map, 
or 40% fewer than the 44 sources in the original map. 
Therefore a number of sources do become confused at 
lower resolution. The basic source properties for the orig- 
inal and degraded-resolution samples, including angu- 
lar deconvolved sizes, axis ratios, mass distribution, and 
mean densities, are compared in Figure [TO) Here the an- 
gular deconvolved size is defined as 9dec = sj S^eas " -^^i 
where 9meas is the geometric mean of the measured mi- 
nor and major axis FWHM sizes and B is the pointing- 
smeared beam FWHM (32.5" and 62" for the original 
and degraded resolution maps, respectively). 

We note that the average source size in the degraded- 
resolution map is nearly twice that in the original Ophi- 
uchus map (average angular deconvolved size of 61" ver- 
sus 98"). Such a large size difference cannot be fully 
accounted for by blending of sources, as even isolated 
sources show the same effect. This behavior provides 
clues to the intrinsic intensity profile of the sources. For 
example, a gaussian intensity profile or a solid disk of 
constant intensity will both have measured deconvolved 
sizes that are similar in maps with 31" and 62" beams. 
Conversely, a power law intensity profile will have a larger 
measured size in the 62" resolution map. 

The ratio of angular deconvolved size to beam size 
{0dec/9mb\ Figure [TT]) arc similar for the degraded- 
resolution (median 0dec/0mb = 1-7) and original (me- 
dian 0dec/Omb = 1-5) samples, further evidence for power 
law intensity profiles. An intrinsic gaussian or solid 
disk intensity profile will result in Odec/^mb values in 
the degraded-resolution map that are approximately half 
those in the original map, while a 1/r^ intensity profile 
results in similar values in the degraded-resolution and 



original map (0.9 versus 1.3). We discuss source profiles 
further in Section [5T3l below. 

Sources in the degraded-resolution map appear slightly 
more elongated, with an average axis ratio at the half- 
maximum contour of 1.3 ± 0.2 compared to 1.2 ± 0.2 for 
the original map. Larger axis ratios are expected for 
blended sources in a lower-resolution map. The slope of 
the mass distribution is not significantly changed: 2.0 ± 
0.4 for the degraded-resolution sample compared to 2.1 ± 
0.3 for the original sample. The factor of two increase in 
deconvolved sizes do create lower mean densities in the 
degraded-resolution sample, however (median density of 
1.6x 10^ cm~^ compared to 5.8x 10^ cm~'^ for the original 
sample). The effect of resolution on mean densities will 
be discussed further in § 15.41 

Given the small number of sources in the degraded- 
resolution map, we carry out the three cloud comparison 
below using the original Ophiuchus map to mitigate un- 
certainties from small number statistics. Any notable 
differences between the original and degraded-resolution 
Ophiuchus results will be discussed where appropriate. 

5.3. Physical Implications of Source Sizes and Shapes 

To compare sources in clouds at different distances, 
we first look at the linear deconvolved size: Ddec = 
d^J Orneas ~ whcrc d is the cloud distance, Omeas is 
measured angular size, and B is the beam size. The dis- 
tributions of linear deconvolved source sizes for the three 
clouds are shown in Figure fT2l left. Sources in Ophiuchus 
have smaller deconvolved sizes than those in Perseus or 
Serpens by almost a factor of two, with mean vales of 
7.5 X 10^ AU in Ophiuchus compared to 1.2 x lO'' AU 
in Serpens and 1.5 x lO"' AU in Perseus. There is a sys- 
tematic uncertainty in the deconvolved size associated 
with the uncertainty in the effective beam size, which 
becomes larger with larger pointing errors. The effective 
beam in any of the three clouds may be as large as 35", 
which would decrease deconvolved sizes negligibly, by up 
to 10^ AU depending on the distance and measured size. 

While there are possible physical explanations for in- 
trinsic size differences, for instance a denser medium with 
a shorter Jeans length should produce smaller cores on 
average, we are more likely seeing a consequence of the 
higher linear resolution in Ophiuchus, as discussed in 
§ 15.21 Thus cores in Serpens and Perseus would likely 
appear smaller if observed at higher resolution, and mea- 
sured linear deconvolved sizes should be regarded as up- 
per limits. To reduce the effects of distance, we exam- 
ine the ratio of angular deconvolved size to beam size 
{Odec/Omb\ Figure [HI right). We found in § [O] that 
Odec/Qmb does not depend strongly on the linear resolu- 
tion, but does depend on the intrinsic source intensity 
profile. 

If the millimeter sources follow power law density dis- 
tributions, w hich do not have a well defined size, then 
lYoung et al.l ([2003, hereafter Y03) show that 9dec/9mb 
depends on the index of the power-law, and not on the 
distance of the source. So, for example, if sources in 
Perseus and Ophiuchus have the same intrinsic power 
law profile, the mean Odec/dmb should be similar in the 
two clouds, and the mean linear deconvolved size should 
be twice as small in Ophiuchus because it lies at half the 
distance. This is precisely the behavior we observe, sug- 
gesting that many of the detected 1.1 mm sources have 
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power law density profiles. 

Considering that a number of the 1.1 mm sources have 
internal luminosity sources (M. Enoch et al. 2007, in 
preparation), and that protostellar envelopes are ofte n 
well described by power law profiles (Shirley et al. 2002, 
Y03), this is certainly a plausible scenario. According 
to the correlation between Odec and density power law 
exponent p found by Y03, median 9dec/0mb values of 1.7 
in Perseus, 1.5 in Ophiuchus and 1.3 in Serpens would 
imply average indices of p = 1.4, 1.5, and 1.6 respec- 
tively. These numbers are consistent with mean p values 
found from radiative trans fer modeling of Cl ass and 
Class I envelopes {p ~ 1.6. IShirlev erall[200l Y03), al- 
though the median for those samples is somewhat higher 
(p ~ 1.8). Note that source profiles could deviate from a 
power-law on scales much smaller than the beam size, or 
on scales larger than our size sensitivity (200") without 
affecting our conclusions. 

Perseus displays the widest dispersion of angular sizes, 
ranging continuously from 1 — 4,9 mb- By contrast, more 
than half the sources in Serpens and Ophiuchus are 
within O.50mh of their respective mean values. Although 
there is a group of Ophiuchus sources at large sizes in 
Figure [121 note that the degraded-resolution Ophiuchus 
sample displays a very narrow range of sizes (Figure fTTj). 
similar to Serpens. The observed size range in Perseus 
would correspond to a wide range of power law indices, 
from very shallow (p ^ 1) to that of a singular isother- 
mal sphere [p = 2). A more likely possibility, however, 
is that sources with large 9 dec 1 9 mh do not follow power 
law density profiles. 

The axis ratio at the half-maximum contour is a simple 
measure of source shape. Figure [T51 shows the distribu- 
tion of 1.1 mm source axis ratios in the three clouds. Our 
simulations suggest that axis ratios of up to 1.2 can be in- 
troduced by the data reduction (Paper I) , so we consider 
sources with an axis ratio < 1.2 to be round, and those 
with a ratio > 1.2 to be elongated. Sources in Ophiuchus 
tend to be round, with a mean axis ratio of 1.2, but note 
that the mean axis ratio in the degraded-resolution Ophi- 
uchus sample is 1.3. The average axis ratio in Serpens 
is 1.3, and Perseus sources exhibit the largest axis ratios 
with a mean of 1.4 and a tail out to 2.7. 

We found in Paper II that Ophiuchus sources were 
more elongated at the 4cr contour than at the half-max 
contour, as would be the case for round cores embedded 
in more elongated filaments. A similar situation is seen 
in Serpens; the average axis ratio at the 4(T contour is 
1.4 in all three clouds. Thus cores in Perseus are some- 
what elongated on average, while objects in Serpens and 
Ophiuchus appear more round at the half-max contour 
but elongated at the 4cr contour, suggesting round cores 
embedded in filamentary structures. 

In addition to angular sizes, Y03 also note a relation- 
ship between axis ratio and density power law exponent, 
finding that asphcrical sources are best modeled with 
shallower density profiles. The inverse proportionality 
between p and axis ratio demonstrated in Figure 25 of 
Y03 suggests power law indices in all three clouds be- 
tween 1.5 and 1.7. These values are consistent with those 
inferred from the average angular deconvolved source 
sizes, and the wider variation of axis ratios in Perseus 
again points to a larger range in p for that cloud. 



5.4. Densities and the Mass versus Size Distribution 

Mean densities calculated using linear deconvolved 
FWHM sizes appear to be significantly higher in Ophi- 
uchus, where the median of the mean densities of the 
sample is 5.4 x 10^ cm~'^, than in Serpens (median 
density 2.2 x 10^ cm~^) or Perseus (median density 
1.6x 10^ cm^'^), as seen in FigurefTil left. There is a large 
scatter with standard deviation of order twice the mean 
value in all three clouds. Sources in Ophiuchus tend to 
be less massive than in the other two clouds, so the larger 
mean densities can be entirely attributed to smaller de- 
convolved sizes in the Ophiuchus sample, which are sen- 
sitive to the shape of the intrinsic density distribution 
(Y03). As noted in 15.21 linear resolution has a strong 
systematic effect on deconvolved sizes, and consequently 
on mean densities. The median density of the degraded- 
resolution Ophiuchus sample is 1.6 x 10^ cm~'^, similar 
to both Perseus and Serpens. 

We additionally calculate mean densities using the full- 
width at 4(7 size rather than the FWHM size (Figure [Ml 
right) to test the hypothesis that mean density differ- 
ences are largely an effect of how source sizes are mea- 
sured. Using this definition, differences between the 
clouds are less pronounced, with median densities of 
1.2 X 10^ cm~'^ in Perseus, 1.3 x 10^ cm~'^ in Serpens, 
and 2.0 x 10^ cm"'^ in Ophiuchus. These numbers sug- 
gest that source mean densities are less dependent on 
cloud distance when measured at the radius where the 
source merges into the background, rather than at the 
half-max. 

Figure [15] again displays the source total mass ver- 
sus size distribution, using the angular deconvolved size 
{9 dec, in units of the beam size) rather than the linear de- 
convolved size. With the exception of two low mass com- 
pact sources, 1.1 mm sources in Serpens exhibit a wide 
range of masses and a narrow range of sizes. Sources in 
Perseus, in contrast, demonstrate a wide range in both 
mass and size. This difference likely reflects a wider vari- 
ety of physical conditions in the Perseus cloud, including 
the existence of more sources outside the main cluster re- 
gions. In such lower density regions, sources may be more 
extended than in densely populated groups. This idea is 
supported by our observations: Perseus sources within 
the NGC 1333 region are smaller on average (mean size 
of 59± 13") than those in the rest of the cloud (mean size 
of 71 ± 21"). Sources in Ophiuchus have a more bimodal 
distribution, with the majority occupying a narrow range 
of sizes, and a smaller group with 9dec > 59mb- 

5.5. Fragmentation and the Core Mass Distribution 

Differential (dN/dM) core mass distributions (CMDs) 
of 1.1 mm sources in Serpens, Perseus, and Ophiuchus 
are shown in Figure [TBI with those of Perseus and Ophi- 
uchus scaled down for clarity. Error bars reflect la ^/N 
statistical uncertainties. Dashed lines indicate empirical 
50% completeness limits for average sized sources in each 
cloud, which are determined as described in § 15.11 and in 
Paper I. Mass distributions include all 1.1 mm cores in 
each cloud, including those that may be associated with 
embedded protostellar sources. 

The shape of the Ophiuchus and Perseus mass distri- 
butions are quite similar above their respective complete- 
ness fimits {M > 0.5Mq in Ophiuchus and M > O.SM© 
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in Perseus). Fitting a power law [dN/dM (x M^") to 
the CMDs, we find that they both have a best fit slope 
of a = 2.1, although the error is larger on the slope 
for Ophiuchus (ctq = 0.3) than for Perseus (ctq, = 0.1). 
The slope of the Serpens CMD is marginally different 
(a = 1.6±0.2), being flatter than in the other two clouds 
by approximately 2(t. 

The two-sided Kolmogorov-Smirnov test indicates a 
high probability (46%) that the Perseus and Ophiuchus 
mass distributions are representative of the same parent 
population. Conversely, the probabilities that the Ser- 
pens core masses are sampled from the same population 
as the Perseus (prob= 12%) or Ophiuchus (prob= 5%) 
masses are much lower. Although it is possible that the 
lower linear resolution in Perseus and Serpens has lead 
to larger masses via blending, the test we conducted to 
increase the effective beam size in the Ophiuchus map 
did not appreciably change the shape of the Ophiuchus 
mass distribution (see Figure [TU)) . 

If the shape of the CMD is a result of the frag- 
mentation process, then the slope of the CMD can be 
compared to models , e.g. of turbulent fragmentation. 
iPadoan fc Nordlundl (|2002[ ) argue that turbulent frag- 
mentation naturally produces a power law with a — 
2.3, consistent with the slopes we measure in Perseus 
and Ophiuchu s (a = 2.1), but not wit h Serp ens (a — 
1.6). Recentlv lBallesteros-Paredes et al.l (|2006L hereafter 
BP06) have questioned that result, finding that the shape 
of the CMD depends strongly on the Mach number of 
the turbulence in their simulations. The BP06 smoothed 
particle hydrodynamics (SPH) simulations show that 
higher Mach numbers result in a larger number of sources 
with lower mass and a steep slope at the high mass end 
(their Figure 5). Conversely, lower Mach numbers favor 
sources with higher mass, resulting in a smaller number 
of low mass sources, more high mass cores, and a shal- 
lower slope at the high mass end. 

Usin g an analytic argument, iPadoan fc Nordltmdl 
(|2002r ) also note a relationship between core masses and 
Mach numbers, predicting that the mass of the largest 
core formed by turbulent fragmentation should be in- 
versely proportional to the square of the Alfvenic Mach 
number on the largest turbulent scale A4\- Given that 
our ability to accurately measure the maximum core 
mass is limited by resolution, small number statistics, 
and cloud distance differences, we focus here on the over- 
aU CMD shapes. 

To compare our observational results to the simulations 
of BP06, we estimate the sonic Mach number Ai = (Jv/cg 
in each cloud. Here a^, is the observed rms velocity dis- 
persion, Cs — \JkT / jimH is the isothermal sound speed, 
and /X — 2.33 is the mean molecular weight per parti- 
cle. Large ^^CO maps of Perseus and Ophiuchus ob- 
served with FCRAO at a resolution of 44" are publicly 
available as part of "The COMP LETE Survey of Star 
Forming Regions"^ (COMPLETE: lGoodman et al.l[200l 
iRjdge et al.l 120061 ). Average observed rms velocity dis- 
persions kindly provided by the COMPLETE team are 
(Tt, = 0.68 km s~^ in Perseus, — 0.44 km s~^ in Ophi- 
uchus, and cr„ = 0.92 km s~^ in Serpens (J. Pineda, per- 
sonal communication) . These were measured by masking 
out all positions in the map that have peak temperatures 

http://cfa-www.harvard.odu/COMPLETE/ 



with a S/N less than 10, fitting a Gaussian profile to each, 
and taking an average of the standard deviations. 

We note that the value of tr^ = 0.68 km s^^ for 
Perseus is smaller than a previous measurement of Uy — 
2.0 km s"i based on AT&T Bell Laborator y 7 m ob- 
serva tions of a similar area of the cloud (Padoan et al.l 
I2003D . The smaller value derived by the COMPLETE 
team is most likely a consequence of the method used: a 
linewidth is calculated at every position and then an av- 
erage of these values is taken. In contrast to calculating 
the width of the averaged spectrum this method removes 
the effects of velocity gradients across the cloud. The dif- 
ferent resolutions of the surveys (44" and 0.07 km s^^ for 
the COMPLETE observations, 100" and 0.27 km s^^ for 
the Padoan et al. observations) may also play a role. 
Sources of uncertainty in the linewidth measurement are 
not insignificant, and include the possibility that ^^CO is 
optically thick to an unknown degree, and the fact that 
line profiles are not necessarily well fit by a gaussian, 
especially in Perseus where the lines sometimes appear 
double-peaked. 

Assuming that the sound speed is similar in all three 
clouds, the relative velocity dispersions suggest that tur- 
bulence is more important in Serpens than in Perseus or 
Ophiuchus, by factors of approximately 1.5 and 2, re- 
spectively. Mach numbers calculated using the observed 
(Ti, and assuming a gas kinetic temperature of 10 K are 
M = 4.9 (Serpens), M = 3.6 (Perseus), and M = 2.3 
(Ophiuchus). We focus on Serpens and Ophiuchus, as 
they are the most different. Serpens is observed to have 
a higher Mach number than Ophiuchus, but the CMD 
indicates a larger number of high mass cores, a shallower 
slope at the high mass end, and a dearth of low mass 
cores compared to Ophiuchus. This result is contrary to 
the trends found by BP06, which would predict a steeper 
slope and more low mass cores in Serpens.^ 

One core of relatively high mass is measured in the 
degraded-resolution map of Ophiuchus, as is expected for 
blending at lower resolution (Figure (TO]), making the dif- 
ference between the Serpens and Ophiuchus CMDs less 
dramatic. The degraded-resolution CMD for Ophiuchus 
still has a steeper slope (a = 2.0 ± 0.4) than the Ser- 
pens CMD, however, and a larger fraction of low mass 
sources: 81% of sources in the degraded-resolution Ophi- 
uchus sample have M < 2Mq compared to 49% in the 
Serpens sample. 

More accurate measurements of the Mach number and 
higher resolution studies of the relative CMD shapes will 
be necessary to fully test the BP06 prediction, given that 
uncertainties are currently too large to draw firm conclu- 
sions. As both numerical simulations and observations 
improve, however, the observed CMD and measurements 
of the Mach number will provide a powerful constraint 
on turbulent star formation simulations. 

Comparing the shape of the CMD to the stellar ini- 
tial mass function (IMF) may give insight into what de- 
termines final stellar masses: the initial fragmentation 

* It should be noted that the shape of the CMDs presented by 
BP06 depend on the type of numerical code used, and not on the 
cloud Mach number alone. Here we have compared to the BP06 
SPH simulations (their Figure 5); for the BP06 total variation di- 
minishing (TVD) method (their Figure 4), a larger number of low 
mass cores are again seen for higher Mach numbers, but the CMD 
slopes at the high mass end are independent of Mach number. 
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into cores, competitive accretion, or feedback pro cesses. 
The shape of the local IMF is still uncertain (jScalo I 
|2005() . but recent work has found evidence for a slope 
of a = 2.5 — 2.8 for stellar masses M > IMq, somewhat 
steeper than the slopes we obs erve in all three clouds. 
For example, Reid et al.l (j2002) find a = 2.5 above 0.6 
Mq, and a = 2.8 above 1 Mq. 
a 2.7 (M > 1 Mq), while 



Chabrieil (l2003f) s uggests 



Schrdder fc Pagel ll2003h 



finds a 2.7 for 1.1 < M < 1.6 Mq and a = 3.1 for 
1.6 < M < 4 Mq. For reference, the Salpeter IMF has a 
slope of a = 2.35 (|Salpctcr 1955). 

For comparison to the IMF, we would ideally like to 
construct a CMD that includes starless cores only, so that 
it is a measure of the mass initially available to form a 
star. Although we cannot separate prestellar cores from 
more evolved objects with millimeter data alone, we are 
currently comparing the Bolocam maps with c2d Spitzer 
Legacy maps of the same regions, which will allow us to 
distinguish prestellar cores from those with internal lu- 
minosity sources (M. Enoch et al. 2007, in preparation). 

5.6. Clustering 
We use the two-point correlation function. 



w[r) = ^ ^ — 1, 



(5) 



Hrir) 

as a quantitative measure of the clustering of cores. Here 
Hs(r) is the number of core pairs with separation be- 
tween log{r) and log{r) ~\- dlog{r), and Hr{r) is similar 
but for a random distribution (see Paper 1). Thus w(r) 
is a measure of excess clustering over a random distribu- 
tion, as a function of separation. 

Figure [17] plots the cloud two-point correlation func- 
tions, with the best fit power laws, w(r) oc r", shown. 
The average linear source FWHM size is indicated, as 
well as the minimum possible separation (the beam size), 
and the linear map size. The best fit slopes to w(r) for 
Serpens {n = -1.5 ± 0.2), Perseus (n = -1.18 ± 0.06), 
and Ophiuchus {n = —1.5 ± 0.1) are consistent within 
2(7, but there is an indication both from the slopes and 
a visual examination of the plot that w(r) falls off more 
steeply in Serpens and Ophiuchus than in Perseus. A 
shallower slope suggests that clustering remains strong 
out to larger scales in Perseus than in the other two 
clouds. 

We do find that a broken power law with slopes of 
n = -0.75 (1.3 X 10"* < r < 6.3 x lO-^AU) and n = 
-3.3 (r > 6.3 X 10"^ AU) is a better fit to the Serpens 
correlation function than a single power law, indicating 
that clustering remains strong to intermediate scales, as 
in Perseus, but then drops off quickly. Also shown in 
Figure [17] is w{r) for the degraded-resolution Ophiuchus 
sample (gray hatched curve). The degraded-resolution 
sample is fit by a shallower slope {n ~ — 1.2±0.3) because 
there are fewer sources at small separations, but it is 
broadly consistent with the original Ophiuchus sample. 

Our general conclusion that clustering of 1.1 mm 
sources remains strong over larger scales in Perseus than 
in Serpens and Ophiuchus can also be reached by visu- 
ally examining each map. Perseus has highly clustered 
regions spread over a larger area as well number 
of more distributed sources, whereas Serpens and Ophi- 
uchus have a single main cluster with fewer small groups 
spread throughout the cloud. 



A few caveats should be noted here. First, the fact that 
clustering seems to extend over larger scales in Perseus 
could simply be due to the fact that Perseus has more 
widely separated regions of star formation, and the phys- 
ical association of all these regions has not been firmly 
established. Second, the map of Serpens covers a much 
smaller linear area than that of Perseus (30 pc^ versus 
140 pc^); we need to be sure that the steepening of the 
slope in Serpens is not caused by the map size. This was 
confirmed by taking a piece of the Perseus map equal 
in size to the Serpens map and recalculating w{r). Al- 
though the value of the slope changed slightly, the smaller 
map size did not cause the slope to steepen at large sep- 
arations. A more serious issue is that the overall ampli- 
tude, but not the slope, of w{r) depends on how large an 
area the random distribution Hr{r) covers. We choose 
each random distribution such that the largest pair sepa- 
ration is similar to the largest pair separation in the real 
data. For Serpens and Ophiuchus this means that the 
random distribution does not cover the entire observed 
area. 

Finally, we look at two other measures of clustering: 
the peak number of cores per square parsec, and the me- 
dian separation of cores. Although the median separa- 
tion of cores is much smaller in Ophiuchus (6.2 x 10"^ AU) 
than in Perseus (2.6 x 10^ AU) or Serpens (2.3 x 10"* AU), 
the median separation in the degraded-resolution Ophi- 
uchus map (2.6 x 10** AU) is consistent with the other two 
clouds. Imposing a uniform flux limit of 75mJy x (250pc-j2 
across all three clouds, equivalent to the 5a limit in the 
shallowest map (Perseus), does not significantly change 
these results. 

We calculate the number of cores within one square 
parsec at each point in the map, and take the peak value 
to be the peak number of cores per square parsec. The 
peak values are 12 pc~^ in Serpens, 22 pc~^ in Perseus, 
and 24 pc~^ in the original Ophiuchus map. Differences 
in linear resolution and completeness both have an ef- 
fect in this case; the peak number in Ophiuchus falls to 
20 pc~^ when using a uniform flux limit, and to 12 pc~^ 
in the degraded-resolution map. The peak number of 
cores per parsec provides further evidence that cluster- 
ing is stronger in Perseus than in the other two clouds. 

5.7. Relationship to Cloud Column Density 

In contrast to the extinction map, which is a measure 
of the general cloud (line-of-sight averaged) column den- 
sity, the 1.1 mm map is sensitive only to regions of high 
volume density (see § 15. ip . A comparison of the two 
tells us, therefore, about the relationship between dense 
star-forming structures and the column density of the 
larger-scale cloud. A visual comparison of the 1.1 mm 
maps of each cloud with visual extinction maps derived 
from the reddening of background stars (e.g. Figure [5]) 
suggests that 1.1 mm cores are generally found in regions 
of the cloud with high Ay ■ 

Figure [18] quantifies the relationship between dense 
cores and the surrounding cloud column density by plot- 
ting the cumulative fraction of 1.1 mm cores in each cloud 
as a function of cloud Ay- In all three clouds the ma- 
jority of cores are found at high cloud column density 
{Av > 7™). The Ay levels above which 75% of cores are 
found in each cloud are indicated by thin lines in Fig- 
ure [TS] 75% of 1.1 mm cores in Perseus, Serpens, and 



1.1 mm Surveys of Serpens, Perseus, and Ophiuchus 



11 



Ophiuchus are found at visual extinctions of Ay > 8™, 
> 15™, and Ay > 23", respectively. Ahhough there 
is not, in general, a strict extinction threshold for finding 
cores, below these Ay levels the likelihood of finding a 
1.1 mm core is very low. Only in Ophiuchus does there 
appear to be a true Ay threshold; only two cores are 
found at Ay < 17™ in this cloud. 

The cumulative distribution for the degraded resolu- 
tion Ophiuchus sample is nearly identical to the origi- 
nal Ophiuchus sample, but shifted to lower Ay by 1- 
2 mag. For the degraded resolution map, 75% of cores 
are found at Ay > 20™. Based on this result and 
the value for the original Ophiuchus sample, we adopt 
Ay 20 - 23™ for the 75% level in Ophiuchus. Re- 
quiring that all cores in the Perseus and Serpens sam- 
ple meet the detection threshold for Ophiuchus (approx- 
imately 110 mJy) changes the distributions negligably. 
Only for Serpens is the 75% level increased slightly, from 
Ay ~ 15™ to 16™. 

Cloud to cloud differences could indicate variations 
in the core formation process with environment, differ- 
ing degrees of sub-structure in the clouds, or varying 
amounts of foreground extinction. Note that in Papers I 
and 11 we found extinction "thresholds" for Perseus and 
Ophiuchus, using a different analysis that looked at the 
probability of finding a 1.1 mm core as a function of Ay, 
of Ay ~ 5™ and Ay ~ 9™, respectively. Those values 
were derived using 2MASS-only Ay maps, rather than 
t he c2d extinction maps used here. 

iJohnstone et al.l ()2004f ) have suggested an extinction 
threshold for forming dense cores in Ophiuchus at Ay ^ 
15™, based on the lowest Ay at which SCUBA cores 
were observed to have sizes and fluxes consistent with 
those of stable Bonnor-Ebert spheres. We find that most 
Bolocam cores in Ophiuchus are found at even higher 
extinctions (75% at Ay > 20 — 23™); the discrepancy 
likely arises frorn differ ences in the extinction maps used. 
IJohnstone et al.l (|2004f) used a 2MASS-derived extinction 
map, while the c2d extinction map used here includes 
IRA C data as well and pr obes somewhat higher Ay val- 
ues. iHatchell et al.l ()2005| ) find no evidence for an extinc- 
tion threshold_Jn_^erseusja^^ / xm SC UBA data, 
but iKirk. Johnstone fc Di Francescol ()2006l ). also using 
SCUBA data, find a threshold of Ay - 5™. Despite 
differences in detail, the relative 75% Ay levels found 
here (highest in Ophiuchus and lowest in Perseus) are 
consistent with the relative values of previous extinction 
threshold measurements in Perseus and Ophiuch us. 

An extinction threshold has been predicted bv iMcKe^ 
(|1989| ) for photoionization-regulated star formation in 
magnetically supported clouds. In this model star for- 
mation is governed by ambipolar diffusion, which de- 
pends on ionization levels in the cloud. Core collapse 
and star formation will occur only in s hielded regions of 
a mol ecular cloud where Ay > 4 — 8™. IJohnstone et al.l 
(|2004f) note that it is not clear how turbulent models 
of star formation would produce an extinction threshold 
for star-forming objects, and interpret such a threshold 
as evidence for magnetic support. 

5.8. Efficiency of Forming Cores 

Another interesting measure of global cloud conditions 
is the fraction of cloud mass contained in dense cores. Ta- 
ble[3]lists the cloud area and cloud mass within increasing 



Ay contours, calculated from the c2d visual extinction 
maps, as well as the total mass in cores within the same 
extinction contour. Cloud masses are calculated from 
the c2d Ay maps 13.1.31 Eq. [2]). The total cloud mass 
within Ay > 2™ is 3470 Mglor Serpens, 7340 Motor 
Perseus, and 3570 Motor Ophiuchus. 

Within a given Ay contour, the mass ratio is defined 
as the ratio of total core mass to total cloud mass, and 
is a measure of the efficiency of core formation at that 
Ay level. For example, the mass ratio at the Ay — 2™ 
contour is equivalent to the fraction of total cloud mass 
that is contained in dense cores. In each cloud 1.1 mm 
cores account for less than 5% of the total cloud mass. 
The mass ratios at Ay — 2™ are similar in all three 
clouds: 3.8% in Perseus, 2.7% in Serpens, and 1.2% in 
Ophiuchus. If we restrict ourselves to Ay > 6™ in each 
cloud, which is reasonable given that the Serpens map 
was only designed to cover Ay > 6™, the mass ratio is 
still low in all three clouds (< 10%), and remains higher 
in Perseus (7%) than in Serpens (4%) or Ophiuchus (2%). 
Low mass ratios are consistent with measurements of the 
overall star-formation efficiency of 1 — 2%, which suggests 
that molecular cloud material is relatively sterile (e.g. 
lEyanslll999l: iLeisawitz et al.l ll989D. 

IJohnstone et al. I (|2004D find a mass ratio of 2.5% 
in a survey of approximately 4 deg^ of Ophiuchus, 
quit e similar to our resu lt. Note, however, that 
the IJohnstone et al.l (|2004( ) core masses (and mass ra- 
tio) should be multiplied by a factor of 1.5 to com- 
pare to our values, due to differences in assumed val- 
ues of Td, Ki,, and d. In a similar analysis of 3.5 
deg^ in Perseus, also using S CUBA 850 ^m data, 
iKirk. Johnstone &: Di Francescol (|2006i) found a mass ra- 
tio for Ay > of only 1%. The difference between this 
result and our value arises primarily from the smaller 
core masses of those authors, which should be multiplied 
by 2.5 to compare to ours. 

In all three clouds the mass ratio rises with increasing 
Ay contour, indicating that in high extinction regions 
a greater percentage of cloud mass has been assembled 
into cores, consistent with the idea that star formation 
is more efficient in dense regions. Although this is an 
intuitively obvious result, it is not a necessary one. If, 
for example, a constant percentage of cloud mass were 
contained in dense cores at all column densities, there 
would be a large number of dense cores lying in a low Ay 
background. On the other hand, a molecular cloud might 
consist of large regions of uniformly high extinction in 
which we would find no 1.1 mm cores because there is no 
sub-structure, and millimeter cores require high density 
contrast (see § 15. ip . 

At higher Ay levels, mass ratios vary considerably from 
cloud to cloud. The mass ratio remains fairly low in 
Ophiuchus, with a maximum value of 9% at Ay = 30™. 
In contrast, the mass ratio rises rapidly in Serpens to 
65% at Ay = 30™, which may suggest that Serpens has 
formed cores more efficiently than Ophiuchus at high Ay . 

6. SUMMARY 

This work completes a three-cloud study of the mil- 
limeter continuum emission in Perseus, Ophiuchus, and 
Serpens. We examine similarities and differences in the 
current star formation activity within the clouds us- 
ing large-scale 1.1 mm continuum maps completed with 
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Bolocani at the CSO. In total, our surveys cover nearly 
20 deg^ with a resolution of 31" (7.5 deg^ in Perseus, 10.8 
deg^ in Ophiuchus, and 1.5 deg^ in Serpens), and we have 
assembled a sample of 200 cores (122 in Perseus, 44 in 
Ophiuchus, and 35 in Serpens). Point mass detection 
limits vary from approximately 0.1 to 0.2 Mq depend- 
ing on the cloud. The results presented here provide an 
unprecedented global picture of star formation in three 
clouds spanning a range of diverse environments. 

These Bolocam 1.1 mm observations naturally select 
dense cores with n > 2 x lO'' cm~'^ and density contrast 
compared to the background cloud of at least 30 — 100. 
We test instrumental biases and the effects of cloud dis- 
tance by degrading the resolution of the Ophiuchus map 
to match the distance of Perseus and Serpens. We find 
that linear resolution strongly biases measured linear de- 
convolved source sizes and mean densities, but not the 
mass distribution slope. Angular deconvolved sizes are 
less strongly affected by cloud distance. 

Rather than a true physical difference, the small mean 
linear deconvolved sizes in Ophiuchus (0.8 x 10** AU) 
compared to Perseus (1.5 x 10** AU) and Serpens (1.2 x 
10'* AU) are likely a result of observing sources with 
power law density profiles, which do not have a well de- 
fined size, at a distance of 125 pc in Ophiuchus versus 
250 pc in Perseus and 260 pc in Serpens. The observed 
mean angular deconvolved sizes and axis ratios in each 
cloud suggest average power law indices ranging from 
p= 1.4 to 1.7 (Y03). 

Sources in Perseus exhibit the largest range in sizes, 
axis ratios, and densities, whereas sources in both Ser- 
pens and Ophiuchus display a fairly narrow range of 
sizes for a large range of masses. We suggest that this 
is indicative of a greater variety of physical conditions 
in Perseus, supported by the fact that Perseus contains 
both dense clusters of millimeter sources and more iso- 
lated distributed objects. A wide range in angular de- 
convolved sizes may also imply a range in the power law 
index of source profiles in Perseus (Y03). 

The slope of the clump mass distribution for both 
Perseus and Ophiuchus is a = 2.1, marginally differ- 
ent than the Serpens slope of a = 1.6. Only Perseus 
and Ophiuchus are consistent within the substantial er- 
rors with the stellar initial mass function (a ^ 2.5 — 2.8) 
and with the slope predicted for turbulent fragmentation 
{a = 2.3) by Padoan & Nordlund (2002). 

Turbulent fragmentation simulations by BP06 predict 
that higher cloud Mach numbers should result in a large 
number of low mass cores, and low Mach numbers in a 
smaller number of higher mass cores. Given the mea- 
sured Mach numbers of = 4.9 in Serpens, 3.6 in 
Perseus and 2.3 in Ophiuchus, our observed core mass 
distribution (CMD) shapes are inconsistent with the tur- 
bulent fragmentation prediction from BP06. We cannot 
rule out a turbulent fragmentation scenario, however, 
due to uncertainties in the observations and in our as- 
sumptions. 

We argue that clustering of 1.1 mm sources remains 
stronger out to larger scales in Perseus, based on the 
slope of the two-point correlation function (-1.5 in Ser- 



pens and Ophiuchus, and -1.2 in Perseus). This result 
is supported by the fact that the peak number of cores 
per square parsec is larger in Perseus (22 pc~^) than in 
Serpens (12 pc~^) or the degraded-resolution Ophiuchus 
map (12 pc^^). 

Finally, we examine relationship between dense cores 
and the local cloud column density, as measured by vi- 
sual extinction {Ay). Extinction thresholds for star for- 
mation have been s uggest ed based on both theory and 
observation (jMcKed fl989l: iJohnstone et al.ll2004f) . Al- 
though in general we do not observe a strict Ay thresh- 
old, dense 1.1 mm cores do tend to be found at high Ay- 
75% of cores in Perseus are found at Ay > 8™, in Ser- 
pens at Ay > 15™, and in Ophiuchus at Ay > 20-23™. 
Our results confirm that forming dense cores in molecu- 
lar clouds is a very inefficient process, with 1.1 mm cores 
accounting for less than 10% of the total cloud mass in 
each cloud. This result is consistent with measurements 
of low star formation efficiencies of a few percent from 
studie s of t he stellar content of molecular clouds (e.g. 
Evansi ll999D . 

While millimeter-wavelength observations can provide 
a wealth of information about the detailed properties 
of star forming cores as well as insight into the large 
scale physical properties of molecular clouds, they do not 
tell a complete story. Detecting and understanding the 
youngest embedded protostars currently forming within 
those cores requires information at mid- to far-infrared 
wavelengths. The Bolocam maps for all three clouds pre- 
sented here are coordinated to cover the same regions as 
the c2d Spitzer Legacy IRAC and MIPS maps of Ser- 
pens, Perseus, and Ophiuchus. Combining millimeter 
and Spitzer data for these clouds will allow us to sepa- 
rate starless cores from cores with embedded luminosity 
sources and to better understand the evolution of cores 
through the early Class and Class I protostellar phases. 
Synthesis and analysis of the combined data set is cur- 
rently underway (M. Enoch et al. 2007, in preparation). 
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TABLE 1 

Sources Found in Serpens 



ID 




RA (2000) 
(h m s) 


; ; 

Dec (2000) 

(° ' ") 


Peak 
(mJy/beam) 




S/N 


other names 


TV TTT^ 

MIPS source.' 


Bolol 


18 28 23.1 


+00 2fi 34 fi 


95 (12) 


5.0 




N 


Bolo2 


18 28 44.0 


+00 53 02 8 


198 (12) 


9.9 




N 


Bolo3 


18 28 45.8 


+00 51 32.4 


227 (12) 


13.7 


IRAS 18262+0050 


Y 


Bolo4 


18 28 47.2 


+00 50 45 1 


145 (10) 


7.1 




N 


Bolo5 


1 8 28 48 3 


+00 14 51 5 


73 ^9") 


5.8 




N 


B0I06 


1 8 28 50 8 


+00 50 28 fi 


115 (10) 


7.3 




N 


Bolo7 


1 8 28 53 


+00 19 03 fi 


119 (11) 


5.4 




Y 


B0I08 


18 28 55.2 


+00 29 28 


617 (13) 


30.8 


IRAS 18263+0027; MMSl (2) 


Y 


Bolo9 


18 28 55.9 


+00 48 30.3 


141 (10) 


9.1 


N 




1 8 28 5fi fi 

-Lcj ^0 iju . yj 


+00 19 10 4 


107 (11) 


5.7 




N 


Bolol 1 


1 8 28 57 3 


+00 48 Ofi 5 


162 (11) 


10.9 




N 


B0I0I2 


1 8 28 58 4 


+00 47 35 7 


172 (11) 


12.0 




N 


Bolol S 


1 8 29 00 2 


+00 30 1 9 8 


239 ^^ 5"! 


5.5 




Y 


B0I0I4 


1 8 29 07 


+00 30 41 5 


1016 ("14") 


59.2 


TRAS 18265+0028- MMS2 Il\ 


Y 


Bolol ^ 


1 8 29 09 fi 

-LCj Wiy . "J 


+00 31 3fi 9 


626 ^^A\ 


33.0 


MMS3 (2) 


Y 


B0I0I6 


1 8 29 1 3 5 


+00 32 1 2 fi 


175 (13) 


5.8 


N 


B0I0I7 


1 8 29 1 fi 4 


+00 18 15 4 


82 (8) 


6.7 


TRAS 18267+0016 


Y 


B0I0I8 


18 29 19 3 


+00 33 29 1 


104 (12) 


6.0 




N 


B0I0I9 


18 29 31 5 


+00 26 49 3 


279 (14) 


14.8 


MMS4 (2) 


N 


Bolo20 


18 29 31 9 


+01 19 00 9 


337 {^ 3") 


18.3 


Y 


Bolo21 


18 29 43.4 


+00 36 25 2 


129 (11) 


7.7 




N 


Bolo22 


18 29 48 8 


+01 16 50 6 


1694 flfil 


80.7 


SMM 9 dV S68N 


Y 


Bolo23 


1 8 29 50 2 


+01 1 5 24 fi 


3010 f20^ 


138 


SMM 1 (1); FIRS 1 


Y 


Bolo24 


18 29 53.8 


+00 36 10.1 


202 (10) 


17.2 


IRAS 18273+0034 


Y 


Bolo25 


18 29 57.4 


+01 13 14.9 


1980 (23) 


62.7 


SMM 4 (1) 


Y 


Bolo26 


18 29 59.2 


+01 14 07.4 


1331 (19) 


43.0 


SMM 3 (1) 


Y 


Bolo27 


18 30 00.3 


+01 10 37.8 


425 (19) 


14.5 


Y 


Bolo28 


18 30 00.7 


+01 12 56.5 


12fi6 (19) 


34.3 


SMM 2 (1) 


Y 


Bolo29 


18 30 01.0 


+01 11 49.1 


979 (19) 


33.5 


SMM 11 (1) 


Y 


BoloSO 


18 30 02.5 


+01 15 24.5 


319 (16) 


11.3 


SMM 8 (1) 


Y 


Bolo31 


18 30 02.8 


+01 08 38.3 


214 (16) 


10.0 


Y 


Bolo32 


18 30 05.7 


+00 39 32.2 


105 (12) 


5.0 




Y 


Bolo33 


18 30 06.4 


+00 42 37.0 


126 (11) 


8.3 




Y 


Bolo34 


18 30 08.2 


+01 13 11.5 


153 (11) 


7.8 




N 


Bolo35 


18 30 14.7 


+01 13 52.fi 


85 (11) 


5.6 




N 



Note. — Numbers in parentheses are \a errors. The peak flux density is the peak value per beam in the lO" 
pixcl"^ unfiltered map (without the optimal fllter applied). The uncertainty in the peak flux density is the local rms 
bcam~^ and does not include an additional 15% systematic uncertainty from calibration uncertainties and residual 
errors after iterative mapping. Other names listed are the most common idcntiflcations from the literature, and 
are not meant to be a complete list. A 1.1 mm sour ce is considered coin cident with a MIP S source if the position is 
within 60^^ of a 24 /xm source fro m the c2d catalog JHarvey et al.|[20Q7D . References — (l) ICasali, Eiroa, &: DuncanI 
{1993); f2) |Diupvik et al.| (|2006l) 



TABLE 2 
Photometry and Core Properties 



ID 


Flux(40") 


Flux(80") 


Flux(120") 


Total Flux 


Mass (lOK) 


Peak Av 


FWHM 


FWHM 


PA 


(n> 


Morphology -"^ 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(A/q) 


(mag) 


(minor,") 


(major,") 


(°) 


cm ^ 




Bolol 


0.137 (0.016) 






0.25 


(0.03) 


0.65 


(0.07) 


8 


52 


(1.0) 


62 


(1.2) 


55 (8) 


l.Ox 10 


extended, elongated, weak 


Bolo2 


0.282 (0.017) 






0.51 


(0.03) 


1.34 


(0.08) 


17 


51 


(0.5) 


63 


(0.6) 


-57 (4) 


2.1x10^ 


multiple , extended , elongated 


Bolo3 


0.324 (0.015) 






0.382 


(0.019) 


1.00 


(0.05) 


19 


43 


(0.4) 


55 


(0.5) 


-77 (3) 


3.2x 10^ 


multiple, extended, round 


Bolo4 


0.226 (0.013) 






0.28 (0.017) 


0.73 


(0.04) 


12 


43 


(0.5) 


55 


(0.6) 


49 (4) 


2.2x 10^ 


multiple , extended , elongated 


Bolo5 


0.094 (0.013) 






0.16 


(0.02) 


0.41 


(0.06) 


6 


50 


(1.2) 


63 


(1.4) 


-68 (8) 


7x lO** 


extended, round, weak 


B0I06 


0.172 (0.013) 






0.210 


(0.016) 


0.55 


(0.04) 


10 


46 


(0.6) 


62 


(0.8) 


-40 (3) 


1.1x10^ 


multiple, extended, elongated 


Bolo7 


0.162 (0.015) 






0.162 


(0.015) 


0.42 


(0.04) 


10 


46 


(0.8) 


57 


(0.9) 


-30 (7) 


l.Ox 10 


multiple , extended , elongated 


B0I08 


0.933 (0.018) 


1.99 (0.04) 




2.63 


(0.04) 


6.87 


(0.12) 


52 


62 


(0.2) 


74 


(0.2) 


68 (1) 


5.2x 10^ 


multiple , extended , round 


Bolo9 








0.169 


(0.010) 


0.44 


(0.03) 


12 


38 


(0.5) 


61 


(0.7) 


-58 (2) 


1.5x 10^ 


multiple , extended , elongated 


BololO 


0.146 (0.015) 






0.146 


(0.015) 


0.38 


(0.04) 


9 


51 


(0.8) 


54 


(0.9) 


11 (25) 


8x 10** 


multiple , extended , elongated 


Bololl 








0.199 


(0.011) 


0.52 


(0.03) 


14 


37 


(0.4) 


61 


(0.7) 


-59 (2) 


1.8x 10^ 


multiple, round. 


B0I0I2 


f\ c\ A ( f\ f\-\ a\ 

0.246 (0.014) 






0.246 


(0.014) 


0.64 


(0.04) 


15 


39 


(0.4) 


58 


(0.6) 


-65 (2) 


2.2x 10 


multiple , extended , elongated 


B0I0I3 


0.36 (0.02) 






0.36 


(0.02) 


0.94 


(0.05) 


20 


52 


(0.5) 


59 


(0.5) 


-78 (6) 


1.7x 10^ 


multiple, extended, round 


B0I0I4 


1.381 (0.019) 






2.16 


(0.03) 


5.65 


(0.09) 


86 


47 


(0.1) 


54 


(0.1) 


-44 (2) 


1.5x 10® 


multiple , extended , elongated 


B0I0I5 


0.856 (0.018) 






0.856 


(0.018) 


2.24 


(0.05) 


53 


44 


(0.2) 


56 


(0.2) 


-27 (1) 


6.5x 10^ 


multiple , extended , round 


B0I0I6 


0.222 (0.018) 






27 


(0 02) 


0.71 


(0.06) 


15 


42 


(0.5) 


62 


(0.8) 


32 (3) 


1.8x 10^ 


multiple, extended, elongated 


B0I0I7 


0.118 (0.011) 


0.24 (0.02) 


0.4 (0.03) 


0.40 


(0.03) 


1.05 


(0.09) 


7 


80 


(1.1) 


97 


(1.3) 


-66 (5) 


3x10** 


extended , elongated 


B0I0I8 


0.144 (0.016) 






0.25 


(0.03) 


0.65 


(0.08) 


9 


48 


(0.8) 


71 


(1.2) 


-41 (3) 


9x 10** 


multiple, extended, elongated, weak 


B0I0I9 


0.41 (0.018) 


0.85 (0.04) 


1.28 (0.06) 


1.28 


(0.06) 


3.34 


(0.14) 


24 


74 


(0.5) 


89 


(0.6) 


-74 (3) 


1.3x 10^ 


extended , elongated 


Bolo20 


0.494 (0.018) 


1.01 (0.04) 


1.57 (0.05) 


1.57 


(0.05) 


4.11 


(0.14) 


29 


87 


(0.5) 


90 


(0.5) 


9 (11) 


1.2x 10^ 


extended, elongated 


Bolo21 


0.177 (0.014) 


0.38 (0.03) 




0.38 


(0.03) 


0.98 


(0.07) 


11 


53 


(0.7) 


64 


(0.8) 


47 (6) 


1.3x 10^ 


multiple, extended, elongated 


Bolo22 


2.32 (0.02) 






3.78 


(0.04) 


9.86 


(0.10) 


144 


47 


(0.1) 


57 


(0.1) 


73 (1) 


2.3x 10® 


multiple, extended, round 


Bolo23 


3.8 (0.03) 


5.91 (0.05) 




5.91 


(0.05) 


15.44 


(0.14) 


256 


46 


(0.1) 


47 


(0.1) 


39 (4) 


6.1x 10® 


multiple , extended , round 


Bolo24 


0.291 (0.014) 






0.40 


(0.02) 


1.04 


(0.05) 


17 


46 


(0.4) 


55 


(0.5) 


-12 (4) 


2.9x 10^ 


multiple , extended , elongated 


Bolo25 


2.75 (0.03) 






3.19 


(0.04) 


8.33 


(0.10) 


168 


44 


(0.1) 


51 


(0.1) 


79 (1) 


2.9x 10® 


multiple , extended , round 


oOlOZD 


i.y t^U.Uoj 






1.90 


(0.03) 


4.95 


(0.07) 


1 1 Q 


48 


(0.1) 


59 


(0.1) 


01 ^^ \ 
-21 (1) 


1 n ^ 1 n6 
l.UX iU 


multiple, extended, round 


Bolo27 


0.65 (0.03) 






1.18 


(0.05) 


3.09 


(0.12) 


36 


47 


(0.3) 


69 


(0.4) 


89 (1) 


4.6x105 


multiple , extended , elongated 


Bolo28 


1.9 (0.03) 






1.90 


(0.03) 


4.96 


(0.06) 


108 


48 


(0.1) 


59 


(0.1) 


-47 (1) 


1.0x10® 


multiple, extended, round 


Bolo29 


1.41 (0.03) 






2.44 


(0.05) 


6.37 


(0.12) 


83 


47 


(0.1) 


59 


(0.2) 


-86 (1) 


1.4x10® 


multiple, extended, round 


BoloSO 


0.43 (0.02) 






0.43 


(0.02) 


1.13 


(0.06) 


27 


48 


(0.4) 


55 


(0.5) 


33 (5) 


2.6x105 


multiple , extended , elongated 


BoloSl 


0.32 (0.02) 


0.65 (0.04) 




0.65 


(0.04) 


1.70 


(0.11) 


18 


45 


(0.5) 


74 


(0.8) 


-65 (2) 


2.4x10" 


multiple, extended, elongated 


Bolo32 


0.151 (0.016) 


0.34 (0.03) 




0.34 


(0.03) 


0.88 


(0.08) 


9 


57 


(0.9) 


74 


(1.2) 


50 (5) 


8x10* 


extended , elongated 


Bolo33 


0.12 (0.015) 






0.128 


(0.019) 


0.34 


(0.05) 


11 


35 


(0.9) 


39 


(1.1) 


25 (17) 


6.9x10® 


round 


Bolo34 


0.181 (0.015) 






0.216 


(0.019) 


0.56 


(0.05) 


13 


40 


(0.6) 


60 


(0.9) 


-70 (3) 


1.7x10" 


multiple , extended , elongated 


Bolo35 








0.06 (0.011) 


0.16 


(0.03) 


7 


24 


(1.1) 


61 


(2.9) 


-69 (4) 


2.4x10" 


multiple, weak 



Note. — Masses arc calculated according to Eq.f^from the total flux density assuming a single dust temperature of To — 10 K and a dust opacity at 1.1mm of Ki.imni — 0.0114 cm^g"~^. 
Peak A^/ is calculated from the peak flux density as in Eq. [3] FWHM and PAs are from an elliptical gaussian flt; the PA of the major axis is measured in degrees east of north, (n) is 
the mean particle density as calculated from the total mass and the deconvolved average FWHM size. Numbers in parentheses are \a uncertainties. Uncertainties for masses are from 
photometry only, and do not include uncertainties from k, To, or which can be up to a factor of a few or more. Uncertainties for the FWHM and PA are formal fltting errors from the 
elliptical gaussian flt; additional uncertainties of 10 — 15% apply to the FWHM, and ~ 5° to the PA (determined from simulations). 

^ The morphology keywords given indicates whether the source is multiple (within 3' of another source), extended (major axis FW at la > l'), elongated (axis ratio at 4cr > 1.2), round 
(axis ratio at 4fT < 1.2), or weak (peak flux density less than 5 times the rms pixel""'^ in the unfiltered map). 
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TABLE 3 

Cumulative Mass as a Function of Extinction 



Mill. Ay 


Area 


Cloud Mass 


Percent Core Mass 


Percent 


Mass Ratio-'^ 


mag 


(%) 


(Mo) 


(%) 


(Mq) 


(%) 


(%) 




Serpens 


2 


100 


3470 


100 


92 


100 


2.7 


4 


97 


3420 


99 


92 


100 


2.7 


6 


60 


2560 


74 


92 


100 


3.6 


8 


28 


1530 


44 


91 


98 


5.9 


10 


15 


1010 


29 


91 


98 


9.0 


12 


9.4 


740 


21 


89 


96 


12 


14 


6.0 


540 


15 


89 


96 


16 


16 


3.9 


390 


11 


86 


93 


22 


18 


2.6 


290 


8.3 


83 


90 


29 


20 


2.0 


230 


6.7 


79 


85 


34 


25 


1.0 


130 


3.7 


61 


66 


47 


30 


0.4 


60 


1.7 


39 


42 


65 


Perseus 


2 


100 


7340 


100 


278 


97 


3.8 


4 


55 


5460 


74 


273 


96 


5.0 


6 


27 


3620 


49 


258 


91 


7.1 


8 


14 


2370 


32 


244 


86 


10 


10 


8.0 


1620 


22 


205 


72 


13 


12 


4.8 


1140 


15 


161 


57 


14 


14 


3.2 


840 


12 


127 


45 


15 


16 


2.0 


610 


8.3 


112 


39 


18 


18 


1.3 


440 


6.0 


96 


34 


22 


20 


0.9 


340 


4.6 


88 


31 


26 


25 


0.4 


180 


2.5 


53 


19 


29 


30 


0.2 


100 


1.4 


6.1 


2.2 


6.1 


Ophiuchus 


2 


100 


3570 


100 


44 


100 


1.2 


4 


55 


2710 


76 


44 


100 


1.6 


6 


25 


1830 


51 


40 


92 


2.2 


8 


14 


1330 


37 


40 


92 


3.0 


10 


9.0 


1070 


30 


40 


92 


3.7 


12 


7.0 


940 


26 


40 


92 


4.3 


14 


5.6 


820 


23 


40 


92 


4.9 


16 


4.6 


730 


20 


40 


92 


5.5 


18 


3.8 


640 


18 


39 


90 


6.1 


20 


3.1 


560 


16 


37 


84 


6.6 


25 


2.0 


400 


11 


32 


73 


8.0 


30 


1.2 


260 


7.3 


23 


52 


8.8 



Note. — Cloud areas and masses are calculated using the c2d extinction maps (sec 
§ [sj of each cloud and Eq. [2] Core masses are taken from Table [2] 

^ The Mass Ratio is computed from the ratio of total core mass to cloud mass within 
the same contour of Ay^ . 
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Fig. 1. — Bolocam 1.1 mm (thick line) and Spitzer c2d IRAC (thin line) and MIPS (gray line) coverage of Serpens overlaid on the 
ICambresvl l fl999l visual extinction map. The area observed with IRAC was chosen to cover Ay > 6"^ in this portion of the cloud. Our 
Bolocam survey covers the same area as the IRAC and a slightly smaller area than the MIPS observations. 
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Fig. 2. — Bolocam 1.1 mm map of 1.5 deg'^ (31 pc'^ at d = 260 pc) in the Serpens molecular cloud. Bolocam has a resolution of 31", and 
the map is binned to 10" pixel"^. The average Icr rms noise is 9.5 mjy beam~^. 
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Fig. 3. - - Map of the local la rms noise per beam in the Serpens map. Positions of the 35 sources are indicated by white plus symbols. 
The average rms noise is 9.5 mjy beam~^, varying by 18% across the map. The noise is higher near bright sources due to sky subtraction 
residuals. 
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Fig. 4. — Bolocam 1.1 mm map of Serpens with the positions of the 35 sources detected above 5cr indicated by red circles. Inset maps 
magnify the most densely populated source regions, including the well known northern Cluster A, Cluster B to the south, and an elongated 
filament reminiscent of the Bl ridge in Perseus (Paper I). Despite the low rms nois e level reach ed (9.5 mjy beam~^), few sources are seen 
outside the cluster regions. Many bright 1.1 mm sources are associated with YSOs IHarvev et al..,200&) . and all are coincident with 160 /im 
emission. 
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Fig. 5. — Visual extinction (Ay) contours calculated from c2d Spitzer maps, as described in ij 13.1.31 overlaid on the gray-scale 1.1 mm 
map. Contours are Ay = 6,9,12,15,20,25™ and are smoothed to an effective resolution of 2.'5. All bright 1.1 mm cores are found in 
regions of high {> 15™) extinction, but not all high Ay areas are associated with strong millimeter sources. 
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Fig. 6. — Distribution of peak (dashed line) and total (solid line) flux densities of the 35 1.1 mm sources in Serpens. Peak Ay values 
derived from the 1.1 mm peak flux densities using Eq. |3] are shown on the upper axis. The mean peak flux density of the sample is 
0.5 Jy beam~^, the mean peak Av is 40™, and the mean total flux density is 1.0 Jy. The 5(t detection limit of 0.05 Jy (dotted line) is 
relatively uniform across the cloud. 
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Fig. 7. — Distribution of source FWHM minor axis (dashed line) and major axis (solid line) sizes, determined by an elliptical gaussian 
fit. The mean FWHM sizes are 49" (minor axis) and 63" (major axis), and the mean axis ratio (major/minor) is 1.3. 




Fig. 8. — Differential mass distribution of the 35 detected 1.1 mm sources in Serpens for masses calculated with To = 10 K. Dotted 
lines indicate the point source detection limit and the empirically derived 50% completeness limit for sources with the average FWHM size 
of 55". The best fit power law (a = 1.6 ± 0.2) is shown, as well as the best fit lognormal function. 
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Fig. 9. — Completeness as a function of linear deconvolved source size in Serpens, Perseus, and Ophiuclms. Symbols show the distribution 
of source mass versus deconvolved size in each of the three clouds, where the size is the linear deconvolved average FWHM. Curves are 
empirical 50% completeness limits determined from Monte Carlo simulations, and demonstrate the dependence of completeness on source 
size and cloud distance. The beam FWHM of 31" corresponds to approximately 8 X 10'^ AU in Serpens and Perseus, and 4 X 10^ AU in 
Ophiuchus. Error bars for average sized sources near the detection limit in each cloud are also shown, as estimated from the results of 
Monte Carlo simulations and pointing uncertainties of approximately 10". 
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Fk;. 10. — Comparison of the basic source properties for the original (black curve) and degraded-resolution, or convolved, (gray hatched 
curve) Ophiuchus source samples, which contain 44 and 26 cores respectively, (a) Measured angular deconvolved sizes arc larger in the 
degraded-resolution map than in the original map by approximately a factor of two (respective mean values of 98" and 61"). (6) Sources 
tend to be slightly more elongated in the degraded-resolution map, with an average axis ratio of 1.3 it 0.2 compared to 1.2 ± 0.2 in the 
original map. (c) The slope of the CMD is not significantly changed for the degraded-resolution sample, but a number of low mass cores 
are blended into a few higher mass sources, (d) Larger deconvolved sizes lead to lower mean densities for the degraded-resolution sample 
(median values of 1.6 x 10^ cm"^ and 5.8 x 10^ cm"^). 
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Fig. 11. — Ratio of angular deconvolved size to beam size {Sd(ic/Qmb) for the original and degraded-resolution, or convolved, Ophiuchus 
maps. Note that is 31" for the original map and 62" for the degraded-resolution map. Measured Odec/^mb values are similar for the 
degraded-resolution (median Odec/^mh = 1-7) and original (median Bd^cl^mh = 1-5) samples, providing evidence for power law intensity 
profiles. An intrinsic gaussian or solid disk intensity profile will result in Qde,cl^mb values in the degraded-resolution map that are nearly 
half those in the original map, while a l/r'^ intensity profile results in similar values in the degraded-resolution and original map (0.9 vs 
1.3). 
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Fig. 12. — Left: Comparison of the distribution of deconvolved linear source sizes in Serpens (solid line), Perseus (shaded), and Ophiuchus 
(dotted line). Histograms are plotted as the fraction of total sources in the cloud as a function of deconvolved size in AU. Estimates of the 
minimum resolvable source size (dashed lines) are based on smearing of the beam by pointing errors. Although measured angular source 
sizes in Ophiuchus are similar to those in Perseus and Serpens, the deconvolved sizes are much smaller due to the different cloud distances. 
Right: Similar, but sizes are measured as the angular deconvolved size in units of the beam FWHM. For isolate d sources with pow er law 
density profiles, Odec/^mb is inversely proportional to the power law index and is independent of cloud distance IjYoung et al.ll2003IV 
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Fig. 13. — Comparison of the distribution of axis ratios, where the ratio is calculated at the half-max contour. Axis ratios < 1.2 are 
considered round, and > 1/2 elongated, based on Monte Carlo simulations. Sources are primarily round in Ophiuchus and Serpens, with 
mean axis ratio of 1.2 and 1.3, respectively. Perseus exhibits the most elongated sources, with a mean axis ratio of 1.4 and a distribution 
tail extending up to 2.7. 
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Fig. 14. — Comparison of the distributions of core mean densities. Left: The mean density is calculated from the total mass (Table |2j 
and the linear deconvolved source size, assuming spherical cores. Despite the lower masses of cores in Ophiuch us c ompared to Serpens and 
Perseus, mean densities in Ophiuchus are higher on average due to small linear deconvolved sizes (see Figure [T2l . The mean densities in 
the degraded-resolution Ophiuchus sample are more consistent with Perseus and Serpens, suggesting cloud distance has a strong effect. 
Right: Similar, except that the mean density is calculated using the full-width at 4(t size rather than of FWHM size. Differences between 
the clouds are diminished; thus the mean density is less dependent on cloud distance when the source size is measured at the radius where 
the source merges into the background, rather than at the half-max. 
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Fig. 15. — Source total mass versus angular deconvolved size (in units of the beam FWH M) f or Serpens, Perseus, and Ophiuchus. 
The angular size is less dependent on cloud distance than the linear deconvolved source size (§ 15. 2p . Sources in Perseus fill more of the 
parameter space than those in Serpens or Ophiuchus, suggesting a wider range of physical conditions in that cloud. Empirically derived 
50% completeness limits are also shown (lines). 
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Fig. 16. — Comparison of the differential core mass distributions (CMDs) of sources in Serpens, Perseus, and Ophiuchus. The Ophiuchus 
curve has been scaled by 1/10, and the Perseus curve by 1/15 for clarity. Uncertainties reflect \/~N counting statistics only. Vertical 
dotted lines show the 50% completeness limits for average sized sources in each cloud (55" FWHM in Serpens, 68" in Perseus and 69" 
in Ophiuchus). The shape of the Perseus and Ophiuchus distributions are quite similar (a = 2.1) in the region where both CMDs are 
complete. The shape of the Serpens CMD is marginally different (a = 1.6) from the other two clouds. 
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Fig. 17. — Comparison of the correlation function w{r) for the three clouds, with power law fits. For each cloud the average source 
FWHM size, resolution limit, and linear map size in AU are indicated. The best fit s lope to w(r) oc r " is n = —1.5 it 0.2 in Serpens, 
n = —1.18 it 0.06 in Perseus, and n = —1.5 it 0.1 in Ophiuchus. The slope found by John stone et al.l (120001) for a smaller map of Ophiuchus 
is shallower (—0.75). The slope of w{r) derived for the degraded-resolution Ophiuchus sample (hatched gray curve) is also shallower 
(—1.2 it 0.3), but the curve appears consistent with the original sample. Differences in the slope of the correlation function likely trace 
differences in clustering properties of millimeter cores in the clouds. 
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Fig. 18. — Cumulative fraction of 1.1 mm cores as a function of cloud Ay for sources in Perseus, Serpens, and Ophiuchus. The cloud 
Ay is derived from the reddening of background stars using 2MASS and IRAC data as described in § |3] Thin solid lines denote the Ay 
levels above which 75% of dense cores are found in each cloud. More than 75% of 1.1 mm cores lie above Ay ~ 8™ in Perseus, Ay ~ 15™ 
in Serpens, and Ay ~ 23™ in Ophiuchus. Results for the degraded resolution Ophiuchus map are shown as well (thin dotted hne), but do 
not vary substantially from the original Ophiuchus sample, with 75% of cores above Ay 20™ . 



